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Abstract 

o ' 

Q . The large value of the top quark mass implies that the rare top decays 

O ; t^ bW+Z, sW+Z and dW+Z, and t cW+W- and uW+W- , are kine- 

^ I matically allowed so long as rrit > my/ + rnz + ruct^ ~ 171.5 GeV + m^. 

PsJ ' or mt > 2mw + fnu^c ~ 160.6 GeV + ruu^c, respectively. The partial decay 

widths for these decay modes are calculated in the standard model. The par- 
tial widths depend sensitively on the precise value of the top quark mass. The 
\ branching ratio for t — > bW^Z is as much as 1 x 10~^ for mt = 200 GeV, 

and could be observable at LHC. The rare decay modes t cW^W and 
' uW~^W~ are highly GIM-suppressed, and thus provide a means for testing 

■ the GIM mechanism for three generations of quarks in the u, c, t sector. 

Now that the top quark mass is known to be quite large, it is possible to examine the 

Oh! 

^ • question of which rare decay modes of the top are kinematically allowed processes. The 
current CDF and DO average value of the top quark mass = 175 ± 8 GeV [|l],H] implies 
that the decays t bW^Z, sW^Z and dW^Z, are allowed decay modes of the top so long 
as rrit > m^r + mz + 171.5 GeV + nid.. The rare decays t cW^W~ and uW^W~ 

also are allowed if nit 2mw + m^^u ~ 160.6 GeV + m^u. For the present central value of 
the measured top quark mass, all of these processes are occurring at or near threshold, and 
are highly phase space suppressed. The decays t — > cW^W~ and uW^W~ are also highly 
GIM-suppressed, and thus are not likely to be seen at standard model rates. The decays 
t — > diW^ Z, however, are not GIM-suppressed and are potentially observable at the LHC. 
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The partial decay widths for these rare decay modes rapidly increase for larger values of 
the top quark mass, and thus are very sensitive to the precise value of the top quark mass. 
Since the decay widths are proportional to | VJ^. i = 1,2, 3, the rare decay t bW^Z (with 
\Vtb\'^ ~ 1) will dominate unless the value of the top quark is below or nearly at threshold 
for this process. 

We begin with the calculation of the partial decay width r{t hW^ Z) in the standard 
model. The branching ratio for this decay process has been computed previously by Decker, 
Nowakowski and Pilaftsis ^ and Mahlon and Parke 0].* The authors of Ref. included 
the finite widths of the W and Z in their calculation, and found a significant enhancement 
in the decay width near threshold due to finite width effects. There is some disagreement in 
the numerical results of Ref. [Q] and . The numerical results presented here are basically 
consistent with the published results of Ref. 0] in the narrow width approximation. There 
is some numerical difference with Ref. which probably stems from the inclusion of finite 
width effects in that calculation. In addition, explicit analytic formulae for the squared 
amplitude of t — > hW^Z are presented in this work. These formulae do not appear elsewhere 
in the literature, and are useful for more detailed studies of the decay mode. Finally, the 
decay widths for the other rare decay modes t cW^W~ and uW^W~ also are computed. 
A search for these decay modes directly tests CKM unitarity in the w-quark sector. 

I. T ^ BW+Z 

The rare decay t bW^Z proceeds via the three tree-level graphs drawn in Fig. 1. The 
amplitudes for these Feynman diagrams are 




) 



I 



{gt^YPL + gt^l^PR} u{pt), (1) 



*The decay process Q qWZ also was considered in Ref. for very heavy fourth generation 



quarks and exotics with mass > 240 GeV. 
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\f2-mbj 



<Pt), (2) 



As = V,, j {^g cos 9w) 6(V4 (^T^T;^] [d'' ' ^) ^(P^)7^PL«(Pt) x (3) 
[-9^' {ks + PwT + g''' {Pz + hr + g^'' i-pz + PwY] , 
where the four momenta ki, and k^ are given by 

ki ^ Pt- Pz ^ Pb+Pw, 

k2=Pt-Pw =Pb + Pz, (4) 
ks^Pt-Pb^Pw+Pz, 
and the couphngs of the Z boson to the left- and right-handed top and bottom quarks are 

9tL = 2 ~ 3 ^ 



9tn = (^—sin^ew^ 



2 
3 

1 1 

1 



9bi^ = [-- + -sin^ew) , (5) 



9bn = sin^ Ow 

In the above amphtudes, Pl^r stand for the left- and right-handed projectors Pl,r = (1 ^ 
75)/2, and 9w is the weak mixing angle. The amplitude A3 depends on the triple gauge 
vertex W^W~Z. This amplitude has been written in unitary gauge, where there is a 
contribution to the W gauge boson propagator proportional to k^k^/rri^. The amplitude 
also can be written in t'Hooft-Feynman gauge = 1), where this contribution is replaced 
by the exchange of the would-be Goldstone boson of the W. 

The total amplitude is given by ^ = .4i -|- .42 + ^3, and the amplitude squared is 

\A\^ = |A|' + I + I + 2AiA; + 2AiA; + 2A2AI (6) 

where the identities 

A.iA-2 — A.'^A.'^-i 

AiAl = A3AI (7) 
A2AI = AzAl, 



have been used. 

The square amphtude \Ai\'^ is 

.4 \ 



lAr = \vti>\' 



9 



2 cos^ 9w J \ki — ml 



X 



^9. 



tL 



(ki ■ pt) (ki ■ pb) - -kl (pt ■ Pb) 



m| 



m 



{Pt ■ Pz) 
■ {Pb ■ Pw) 



w 



{kl ■ Pb) {kl ■ Pz) - {Pb ■ Pz) 
{kl ■ Pt) {kl ■ Pw) - ^kl {pt ■ Pw) 



{Pt ■ Pz) {Pb ■ Pw) 



W"''Z 



{kl ■ Pz) {kl ■ Pw) - ]jkl {pw ■ Pz) 



-'^'^'>TT't9tL9tR \ {kl • Pb) H 2~ (^1 • Pw) {Pb • Pw) 



m 



w 



{2 2 
{Pt ■Pb) + — {Pt ■ Pz) {Pb • Pz) + {pt ■ Pw) {Pb • Pw) 



+ 



{Pt • Pz) {Pb • Pw) {pw • Pz) 



(8) 



The square amphtude |^2p is related to by the interchanges gt^ ^ 9bL ri ''^f^ ''^bi 
Pt ^ Pb and kl k2, 

\ / 1 \ 2 



\A2\' = \Vt 
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tb 



2cos^9wJ \k2 — ml 



X 



^^9 



bL 



{kl ■ Pb) {k2 ■ Pt) - {Pt ■ Pb) 



mi 



m 



{Pb ■ Pz) 
- {Pt ■ Pw) 



{k2 ■ Pt) {k2 ■ Pz) - 2^2 {Pt ■ Pz) 
{k2 ■ Pb) {k2 ■ Pw) - ]:kl {pb ■ Pw) 



w 



m^m\ 



{Pb ■ Pz) {Pt ■ Pw) 



(/C2 • Pz) {k2 ■ Pw) - ^kl {pw ■ Pz) 



-'^'^nT'l9bL9bR { {k2 • Pt) H 2~ (^2 • Pw) {Pt ■ Pw) 



m 



w 



{2 2 
{Pt •Pb) + — {Pb • Pz) {pt ■ Pz) + {Pb • Pw) {Pt ■ Pw) 
mz my/ 



+ 



m^ml 



{Pb ■ Pz) ipt ■ Pw) {Pw ■ Pz) 



(9) 



The square amphtude \Ai\ is 



\As\'= \V,,\' 



2 J Ui - m 



2 _ ^2 

W, 



X 



4(pt-Pfe)|-3(m^ + m|)+27 ( 



- "^1) - 7^ ("^M/ + ^z) +2(1- 7^) {pw ■ Pz) 



+ 
+■ 



1 2/ 1 



2 \mly 



rrir. 



my m 



w , 



2^ \w? 



w 



rrir 



{pw -pzY 



m 



^ (1-7^) {PwPzf\ 
w^z ^ ' > 



(10) 



+4 [pt ■ Pw) (Pb • Pw^) I -2 - 47 + 272 + 4- (1 + 7) iPw ■ Pz) + (1 + 27 + 72) {p^ ■ pzf \ 

\ T^W fTT'w'^Z J 

+4 ipfPz) ipb-Pz) I -2 + 47 + 272 + A. (1 _ ^) {p^.pz) + -^-^ (1 - 27 + 7^) {pwPzf \ 
+4 [{pt ■ Pw) {Pb ■ Pz) + {Pb ■ Pw) ipt ■ Pz)] I - 6 + 27^ 



+ 



1 



1 



777.7 777? 



27 a 



1 



''Z '"'W / \^z ^w , 

where 7 = (1 — 777|/777^) in unitary gauge. In t'Hooft-Feynman gauge, one obtains the 
same expression with 7 = siv? Oy/ and cos^ replaced by (7772/7771^)^- The terms in |v4.3p 
proportional to 7 are antisymmetric under pf <-> p^, pw Pz and 777^^^ <-> 777^ while the 
terms which are independent of 7 or proportional to 7^ are invariant under this interchange. 
The interference term ^1.42 is 



{pw ■ Pz) 



1 



1 - 



{pw-PzY 
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1 



1 



tb\ 



2cos^ 6w) xkl — ml) ykl — ml 



X 



- 6 {ki ■ k2) {pt ■Pb)-'^ {ki ■ Pt) {k2 ■ Pb) - {ki ■ pb) {k2 ■ Pt) 
{ki ■ /C2) {Pb ■ Pz) {pt ■ Pz) + {ki ■ Pz) {k2 ■ Pz) {pt ■ Pb) + {k2 ■ Pz) {ki ■ Pt) {Pb ■ Pz) 



+ (ki ■ Pz) (pt ■ Pz) (k2 ■Pb)-'^ (Pt ■ Pz) (ki ■ Pb) (k2 ■pz)-'2 (ki ■ pz) ipb ■ Pz) {k2 ■ Pt) 
4 



777 



w L 



{ki ■ k2) (pb ■ Pw) (pt ■ Pw) + (ki ■ Pw) (k2 ■ Pw) (pt ■ Pb) + (^2 • Pw) (ki ■ pb) (pt ■ Pw) 



+ (ki ■ Pw) (Pb ■ Pw) {k2 ■Pt)-'^ (Pb ■ Pw) {ki ■ Pt) {k2 ■ Pw) - 2 (/ci • pw) (pt ■ Pw) {k2 ■ Pb) 
4 



m^m% 



{ki ■ k2) {pw ■ Pz) {- {Pt ■ Pb) {Pw ■ Pz) + {Pb ■ Pz) {Pt ■ Pw) + {Pb ■ Pw) {Pt ■ Pz)} 



+ ipw ■ Pzf {{ki ■ Pt) {k2 ■ Pb) + {ki ■ Pb) {k2 ■ Pt)} 
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+ {Pt • Pb) {pw • Pz) {{h ■ Pz) {h ■ Pw) + {h • Pw) {h ■ Pz)} 
+2 {ki ■ Pz) {k2 ■ Pz) {Pb ■ Pw) {Pt ■ Pw) + 2 (/ci • pw) (/c2 • Pw) {Pb ■ Pz) {pt ■ Pz) 

- {ki ■ Pb) {k2 • Pz) {pt • Pw) {pw ■ Pz) - {k2 ■ Pt) {ki • Pz) {Pb ■ Pw) {pw • Pz) 

- {ki ■ Pt) (fe • Pz) {Pb ■ Pw) {pw ■ Pz) - {k2 ■ Pb) {ki ■ Pz) {pt ■ Pw) {pw ■ Pz) 

- {ki ■ Pt) {k2 ■ Pw) {Pb ■ Pz) {Pw • Pz) - {k2 ■ Pb) {ki • Pw) {Pt • Pz) {Pw • Pz) 

- {ki ■ Pb) {k2 ■ Pw) {pw ■ Pz) {pt ■ Pz) - {k2 ■ Pt) {ki ■ Pw) {Pb ■ Pz) {pw ■ Pz) 

(2 4 

+'^'>T4gti^gbA 3 {ki • p*) H ^ {ki ■ pz) {pt • Pz) r (Pt " Pw) {ki ■ Pw) 

L "^w 

{Pw ■ Pz) [{ki ■ Pw) ipt ■ Pz) + ipt ■ Pw) {ki ■ Pz) - {pw ■ Pz) {ki ■ Pt)] | 



(2 4 

+'^f^htR9bL i 3 (/ca • Pb) H 2" (^2 • Pz) {Pb ■ Pz) 2~ (Pb ' Pw) {k2 ■ Pw) 

I niz "^w 

2 r 

H — 2 — 2 (Pw ■ Pz) {k2 ■ Pw) {Pb ■ Pz) + {Pb ■ Pw) {k2 ■ Pz) - {pw ■ Pz) {k2 ■ Pb) 

ITT'W^Z 

^-2mlmlgt^gb^ \ -5 + ^ {pw ■ Pzf \ J 
I ^w^z ) J 

A1A2 is invariant under the simultaneous interchanges gtj^ r QbL ri ''^t ^ f^bi Pt Pb and 
ki k2- 

The two interference terms AiA\ and A2A\ are 



S'ti { - 4 [(2 - 7) {ki ■ Pw) {pt ■ Pb) + {ki ■ Pb) {Pt ■ Pw) - 3 (/ci • pt) {pb ■ Pw)] 
+4 [(2 + 7) {ki ■ Pz) {pt ■ Pb) + {ki ■ Pt) {Pb -Pz)-^ {ki ■ Pb) {Pt ■ Pz)] 

4 

— 2~ (1 + (^1 ■ Pw) [2 ipt ■ Pw) {Pb ■ Pw) + {Pt ■ Pw) {Pb ■ Pz) 



- {pw • Pz) {pt • Pb) + {Pb • Pw) {pt ■ Pz)] 

4 

+ — (1 - 7) (^1 • Pz) [2 ipt ■ Pz) {Pb ■ Pz) + {pt ■ Pz) {Pb ■ Pw) 
rriz 

- (Pw ■ Pz) (pt ■ Pb) + (Pb ■ Pz) (pt ■ Pw)] 

2 



m 



2 

w 



(3 + 7) {Pw ■ Pz) [{ki ■ Pw) {Pt ■ Pb) - {Pt ■ Pw) {ki ■ Pb) + {Pb ■ Pw) {ki ■ Pt)] 



2 

+ — (3 - 7) {Pw ■ Pz) [{ki ■ Pz) {pt ■ Pb) - {.Pb ■ Pz) {ki ■ Pt) + {pt ■ Pz) {ki ■ pb)] 
rriz 



H — 2 — 2 (1 + 7) {Pw ■ Pz) {Pb ■ Pw) [{ki ■ pw) {Pt ■ Pz) - {Pw ■ Pz) {ki ■ Pt) + {pt ■ Pw) {h ■ Pz)] 
n^w^z 
4 

2 2 (1 ~ 7) {Pw ■ Pz) {Pt ■ Pz) [{h ■ Pw) {Pb ■ Pz) - {Pw ■ Pz) {h ■ Pb) + ipb ■ Pw) {h ■ Pz)] 

fTT-W^Z 

+^rnl9tR I (-4 + 7) (P6 • Pw) + (-2 + 7) ipb ■ Pz) + ^ + ^7) {Pw ■ Pz) {Pb ■ Pw) 

+^ - ^7) {Pw ■ Pz) {Pb ■ Pz) + ^ 2 (1 + 7) {Pw ■ Pzf ipb ■ Pw) \ (12) 



and 



|2 ( 9_ 



^^A-\Vtb\^'-:r][T^--A[T^—^]>< 



m 



4 [(2 - 7) {k2 ■ Pw) {Pt ■ Pb) + {h ■ Pt) {Pb ■ Pw) - 3 (A;2 • Pb) {pt ■ Pw)] 
+4 [(2 + 7) {k2 ■ Pz) {pt ■ Pb) + {k2 ■ Pb) {pt ■Pz)-S {k2 ■ Pt) {Pb ■ Pz)] 

4 

— (1 + 7) {k2 ■ Pw) [2 ipb ■ Pw) {Pt ■ Pw) + {Pb ■ Pw) ipt ■ Pz) 



w 

- (pw ■ Pz) (pt ■ Pb) + (pt ■ Pw) (Pb ■ Pz)] 

4 

— 2 (1 - 7) (^2 • Pz) [2 ipb ■ Pz) {Pt ■ Pz) + {Pb ■ Pz) {Pt ■ Pw) 

- {Pw ■ Pz) (pt ■ Pb) + (pt ■ Pz) (Pb ■ Pw)] 

2 



m 



2 

w 



(3 + 7) {pw ■ Pz) [{k2 ■ Pw) {Pt ■ Pb) - {Pb ■ Pw) {k2 ■ Pt) + {pt ■ Pw) {k2 ■ Pb)] 
(3 - 7) (Pw ■ Pz) [{k2 ■ Pz) {Pt ■ Pb) - {Pt ■ Pz) {k2 ■ Pb) + {Pb ■ Pz) {k2 ■ Pt)] 



4 

2 — 2 (1 + 7) {Pw ■ Pz) {pt ■ Pw) [{k2 ■ Pw) {Pb • Pz) - {pw ■ Pz) {k2 ■ Pb) + ipb ■ Pw) {k2 ■ Pz)] 



4 

2 — 2 (1 ~ 7) {Pw ■ Pz) {Pb ■ Pz) [{k2 ■ Pw) {Pt ■ Pz) - ipw ■ Pz) {k2 ■ Pt) + {Pt ■ Pw) {k2 ■ Pz)] 

I (-4 + 7) {Pt ■ Pw) + (-2 + 7) {Pt ■ Pz) + -\- C- + ^7) {Pw ■ Pz) {Pt ■ Pw) 
y m^r V2 2 / 



1 

+ — 
mi 



~ ^7) {PW ■ Pz) {Pt ■ Pz) + 2^ 2 (1 + 7) {PW ■ Pz)^ {Pt ■ Pw) \ , (13) 



where 7 = 1 — m|/m^ or sin^ 9w- The interference terms ^2-4-3 and ^1^3 are related by 
the interchanges gtj^ a QbL ri ''^t 'f^bi Pt Pbi and ki ^ k2. 

The above square amphtudes have been written in terms of ki and A;2, and pt and p^j 
in order to exhibit the symmetries of the square amphtudes exphcitly. The total square 
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amplitude can be rewritten in terms of the three dot products {pb ■ Pw), {Pb ■ Pz), and 
{Pw ■ Pz), by ehminating pt, ki and k2 in the above formulae. 

The partial width for the decay mode t —>■ bW^Z is given by the three-body phase space 
integral 

T{t^ bW^z) = / dml, dm\, W (M) 

where the invariant square masses mf- = {pi + Pj)^ are defined in terms of the momenta of 
the final particles, and the spin-averaged square amplitude 

W=l\A\^ (15) 

since one averages rather than sums over the top quark spin. 

The partial decay width T {t ^ bW^Z) is plotted in Fig. 2 as a function of the top 
quark mass. The phase space integral was performed numerically for the parameter values 
mw = 80.3 GeV, mz = 91.2 GeV, = 4.5 GeV, sin^ Ow = 0.23 and \Vtb\ = 1. The partial 
width is plotted over the range from = 176 GeV, where the partial width vanishes, 
to rrit = 200 GeV, where the partial decay width is 2.5 x 10^^ GeV. The branching ratio 
BR(t bW^Z) also is plotted as a function of the top quark mass in Fig. 3, assuming that 
the total width of the top quark is dominated by t — > bW^ , 

^ ^ 647r my^ml ^ ' 

\m\ + ml- 2m^ + m'fm'^ + mlm\r - 2m'fml} , (16) 

where 

X{x, y, z) = + + — 2xy — 2yz — 2xz. (17) 

The CKM matrix element \Vtb\^ cancels out of the branching ratio. The branching ratio 
increases from zero for rrit = 176 GeV to 1.0 x 10^^ for rrit = 200 GeV.[| Although this 



^The value of the branching ratio for this value of mt is consistent with the result of Mahlon 
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branching ratio is too small to be observed at the Tevatron, it is large enough to be interesting 
for the LHC which is expected to yield about a million fully-reconstructed top quark events 
per year 0. The observability of this decay mode depends on the precise value of the top 
quark mass. The branching ratio is greater than 10~® for rrit > 187 GeV. The extreme 
sensitivity of the branching ratio to the top quark mass implies that the decay mode could 
be used to extract or bound the top quark mass. The decay mode also is sensitive to the 
presence of the triple gauge vertex W^W~Z with the standard model coupling. 

The rare decay t bW^Z is at threshold for the present central value of the top quark 
mass, so the decays t —>■ sW'^Z and t — > dW^Z could be more important if t — >■ bW^Z is 
kinematically forbidden or just allowed. Alternatively, it might be possible to look at these 
modes by applying a tight cut on the invariant mass of the and Z momenta to exclude 
t bW^Z but not t sW^Z and dW^Z. The partial decay widths for the s and d final 
states can be obtained from the partial decay width for the b mode by replacing mf, by or 
md, and by \Vts\'^ or \Vtd\^. The partial decay width F (t ^ diW^Z) /\Vti\^ is plotted as 
a function of the top quark mass for vanishing m^^- in Fig. 4. The partial width divided by 
the CKM matrix element squared is zero at threshold where = m\Y + mz and increases 
to 2.7 X 10^^ GeV at irtt = 200 GeV. For canonical values of [Visp and |V^dp, these partial 
widths will be too small to be observed at LHG. 

II. T CW+W- 

The rare decay t cW^W~ proceeds through tree-level graphs with intermediate d, s, 
and b quarks, as depicted in Fig. 5. The amplitude for the decay is 

and Parke Q in the narrow width approximation. There is some numerical difference with their 
narrow width results for smaller values of rrit, which probably stems from the inclusion of finite 
width effects proportional to Tw/m]y and Tz/mz in their calculation. 



9 



(18) 



j=d,s,b 

where k = Pc+ Pw- = Pt — Pw+ ■ The spin- averaged amphtude squared is given by 



1 



1 



\V2 J — rrij J — ml^ 



X 



1 



(k-pt) (k-pc) - -k {pc-Pt) 



m 



w 



{PfPw+) 
+ (Pc-Pw-) 



(k-Pc) (k-pw+) - {Pc-Pw+) 



{k ■ Pt) {k ■ pw- 



-k^ iPfPw-) 



(19) 



{PfPw+) {Pc-Pw- 



{k ■ pw+) {k ■ Pw-) - iPw+ ■ Pw-) 



where the factor of 1/2 comes from averaging over the top quark spin. Note that this square 
amphtude can be derived from Eq. (8) or (9). The amphtude squared can be rewritten in 
terms of the three dot products {pc ■ Pw+), {Pc ■ Pw-) a-nd ipw+ ■ Pw-) of the final particle 
momenta by eliminating pt and k in the above formula. 

The partial width for the decay mode t cW^W~ is given by the three-body phase 
space integral 

1 1 



dm23 dmi2 \A\ , 



(20) 



(27r)^ 32m? 

where the invariant square masses m? = (pj + Pj)^ are defined in terms of the momenta of 
the final particles. It is possible to do the m^+p^- integral explicitly, so that the partial 
width is given by an integral over x = = (pc + Pw-)"^ = {Pt — Pw+)'^ = k"^, 



dxi 



1 1 
2 

1 



X 



X — 



^x — m| J 2x 



A2 (a;,m^,m^) A2 (a;,mj^m^) 



X 



(^x + - m^) (. 



(^x — rrif — 



X + ml - rn^ 



) 



(21) 



m 



w 



ml, 



X + — iTiy^.' ] [x — rUf +m 



) 

+ (x - ml - m^) (x + nif - m^) (x — ml + m^) 
10 



m 



w 



(x — ml + m^) {x - ml - m^) {x -m^.^ m^) (x — ml — m^) 



The integrand of Eq. (21) is symmetric under the interchange ml ^ ml, but this symmetry 
is broken by the hmits of integration of the remaining phase space integral. An important 
observation about the decay width is that the width vanishes for = or m| = by CKM 
unitarity, 



j=d,s,b 



(22) 



This GIM suppression can be made manifest by replacing the two d-quark propagators by 



X — m; 



^x — mj J X 



m. 



which implies that the integrand is multiplied by 



2 2 

m^mi 
x'^ 



x{x — m'fj 



(23) 



(24) 



Thus, the final formula for the partial width is 

4 



r{t^ cW+W-) =J^s^3 \ [j=) E ytjy:jV;,V^k / (mj, mi mi, ml, m^) , (25) 
where the integral equals 

7-/2 2 2 2 2\ 

I [mj,mk,m^,mt,mw) 

(x, ml, m^) (x, ml, m^) | + ~ ^^vf) (x + ml - m^) 



r{mt-mw)'^ f 1 
m^mj^ I ax I ^ 



■3/ \X-T^k, 



2x3 



m 



w 



m 



w 



[x — ml — ml^ {x + ml — m^) [x — ml + t^w) 
+ {x - ml - m^) {x + ml - m^) {x -ml + m^) 
(a; - + m^) (x — ml — m^) - + m^) (^x - - 



(26) 



Numerical integration of Eq. (26) (which assumes that GIM suppression is operative) 

shows that the decay width is completely dominated by the contribution with m'j = m\ = 
mf. The partial width is plotted as a function of the top quark mass for m^ — 80.3 GeV, 
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rric — 1.5 GeV, sin^^^ = 0.23 and Vcb — 0.036 — 0.046. The two curves correspond to 
the lower and upper values of the CKM matrix element Vcb- The partial width vanishes at 
threshold where rrit = rric + 2mw, and is at most pa 10~^^ GeV for rrit = 200 GeV. This 
extremely small partial width is a direct consequence of three-family unitarity of the CKM 
matrix in the rt-quark sector. If the GIM suppression condition Eq. (22) is relaxed, the 
integral appearing in Eq. (21) is a factor of 2 x 10^ larger than I {ml, ml, ml, 171^,171^) for 
each value of m| and m^. Thus, it is quite possible that the rare decay t — > cW^W~ occurs 
at an observable level in non-standard model theories. A search for this rare decay mode 
would directly test CKM unitarity of the tc rows of the CKM matrix. 

The partial width for the rare decay t — > vW^W~ can be obtained from the above with 
the replacement c <-> m. The partial width for the up mode is even smaller than for the 
charm mode due to smaller CKM matrix elements. This decay mode can be used to test 
CKM unitarity of the tu rows of the CKM matrix. 

III. CONCLUSIONS 

The partial widths for the rare top decay modes t —>■ bW^Z, sW^Z, dW^Z, cW^W~ 
and uW'^W" have been calculated in the standard model. The decay mode t — > hW^Z 
is potentially observable at LHC rates for top quark masses above 187 GeV, and could be 
used to accurately determine the top quark mass. The decay amplitude also depends on 
the triple decay vertex W^W'Z, and therefore tests for the presence of this coupling and 
its value. The decays t — > cW'^W~ and uW'^W~ are extremely GIM-suppressed in the 
standard model, but may be much larger in non-standard scenarios. A search for these rare 
decay modes tests CKM unitarity of the tc and tu rows of the CKM matrix, 

E VtjVlj-^- (27) 

j=d,s,b 
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FIGURES 

FIG. 1. t ^ hW^Z. Feynman diagrams correspond to the amplitudes Ai-, A2, and ^3. 

FIG. 2. r{t hW+Z) as a function of the top quark mass for mw = 80.3 GeV, mz = 91.2 GeV, 
mb = 4.5 GeV, sin^ 6w = 0.23, and |14fep = 1. The partial decay width vanishes at threshold where 

mt = mb + mw + mz- 

FIG. 3. BR{t — > bW~^Z) as a function of the top quark mass for mw = 80.3 GeV, 
mz = 91.2 GeV, mb = 4.5 GeV, and sin^ 9w = 0.23. The branching ratio vanishes at thresh- 
old where m^ = mb + mw + mz- 

FIG. 4. T{t — s- diW'^Z)/\Vti\'^ as a function of the top quark mass for mw = 80.3 GeV, 
mz = 91.2 GeV, sin^ 9w = 0.23 and m^^ = 0. The partial decay width vanishes at threshold where 
mt = mw + nT-z- This graph is relevant for the decays t — dW^Z and t — sW^Z. 



FIG. 5. t cW+W-. 

FIG. 6. r(t — > cW'^W^) as a function of the top quark mass for mw = 80.3 GeV, 
mc = 1.5 GeV, sin^ 9w = 0.23 and Vcb = 0.036 — 0.046. The partial decay width vanishes at 
threshold where m^ = mc + 2mw- 



15 



Figure 1 
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